%% Transformation Function %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function c1 = trans(c0) % Constrains values of regular coefficients

    c1 = c0;
    
         c1(1) = exp(c0(1))/(1+exp(c0(1))+exp(c0(2)));
         c1(2) = exp(c0(2))/(1+exp(c0(1))+exp(c0(2)));

         c1(3) = exp(c0(3))/(1+exp(c0(3)));
         c1(4) = exp(c0(4))/(1+exp(c0(4)));
         
         c1(6) = -exp(c0(6)); %use Hamilton restriction to identify L-shape regime

         c1(9)  = sqrt(exp(c0(9)));
         c1(10) = sqrt(exp(c0(10)));
   

end % End of transformation function